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A close coupling quantum-mechanical calculation is performed for rotational energy transfer in a 
HD+HD collision at very low energy, down to the ultracold temperatures: T ~ 10~ 8 K. A global 
six-dimensional H2-H2 potential energy surface is adopted from a previous work [Boothroyd et ai, 
J. Chem. Phys., 116, 666 (2002).] State-resolved integral cross sections a^^iy (skin) of different 
quantum-mechanical rotational transitions ij — > i'j' in the HD molecules and corresponding state- 
resolved thermal rate coefficients kij^iji (T) have been computed. Additionally, for comparison, 
H2+H2 calculations for a few selected rotational transitions have also been performed. The hydrogen 
and deuterated hydrogen molecules are treated as rigid rotors in this work. A pronounced isotope 
effect is identified in the cross sections of these collisions at low and ultracold temperatures. 

PACS numbers: 34.50. Cx 34.50.Ez 



I. INTRODUCTION 

The recent creation and investigation of a quantum gas 
of ultracold diatomic molecules [l[ is of great interest in 
many areas of atomic, molecular, optical, and chemical 
physics [H-Q. Research in these fields may have impor- 
tant future applications, for example, in quantum infor- 
mation processing [7H10l| . From a scientific point of view 
the creation of the molecular quantum gas opens new 
doors, for instance, in the experimental and theoretical 
investigation of the cold and ultracold molecular scatter- 
ing and chemical reactions [lll - [l6j . It allows researchers 
to probe the interaction and collisional properties of dif- 
ferent light and heavy molecules in the cold and ultra- 
cold regime: T ~ 10" 4 - 10" 8 K @, S [H EE El • In 
this regime, one can expect many shape resonances in 
the cross sections arising from the van der Waals force 
0, Q. For example, a resonance with a weakly bound 
level near zero collision energy can significantly enhance 
the tunneling effect through a reaction barrier. By align- 
ing and orienting the colliding molecules, the anisotropy 
of the van der Waals forces enables substantial tuning of 
the molecular levels to create such resonances @. 

In this work, the ultracold collision between two 
deuterated hydrogen molecules, i.e. rotational energy 
transfer in HD+HD, is mainly considered. For compar- 
ison similar collision in H2+H2 is also considered. From 
a theoretical point of view the HD+HD system is inter- 
esting because its PES can be derived from the much- 
studied H2+H2 system by adjusting the coordinate of 
the HD-molecule center of mass. Once the symmetry is 
broken in H2-H2 by replacing the H with a D atom in 
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each H 2 we have the precise HD-HD PES. The HD and 
H2 molecules are treated as rigid monomer rotors in this 
work, so we ignore the vibrational degrees of freedom of 
these molecules. 

Because of the small reduced mass and large 
rotational-energy spacing in the HD-HD system, the 
number of states required in the basis set for an accu- 
rate quantum-mechanical calculation should be relatively 
small. The HD+HD system has widely spaced rotational- 
energy levels and, because of the strong anisotropy 
of the intermolecular potential, it has relatively large 
rotational-energy transition probabilities. Since HD is 
a light molecule it can be manipulated easily by an ex- 
ternal electrical field and also, the laser cooling of this 
diatomic molecule seems possible making this system is 
of current experimental interest. 

Surprisingly, such a fundamental and attractive quan- 
tum four-atomic system has not received substantial at- 
tention in previous experimental and theoretical inves- 
tigations. Several molecular-beam studies of HD+HD 
involve the measurement of a few rotational probabilities 
[l9[, integral cross sections for unresolved internal [20j |. 
and rotational energy transfer rates [U [22[ . Neverthe- 
less, there are only a few calculations dealing with the 
rotational excitation in the HD+HD collision, e. g., an 
early modificd-wave-number calculation by Takaya nag i 
[23|, semiclassical calculations by Gelb and Alper |24j . 
and Cacciatore and Billing J25J . 

Hydrogen is otop e effects have often attracted consider- 
able attention |26| . In this work we carry out such consid- 
eration within the HD+HD and H2+H2 systems at high, 
low and ultracold temperatures. In the next section we 
briefly present the quantum-mechanical approach, that is 
used in this work and the PES. In Sec. Ill we present nu- 
merical results for both HD+HD and H2+H2 collisions. 
Additionally, we present a brief discussion of the numer- 
ical convergence of the results. Finally, in Sec. IV we 
present a summary and conclusion. 
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FIG. 1: Four-atomic system (12) + (34) or HD+HD, where 
H is a hydrogen atom and D is deuterium^ represented by 
few-body Jacobi coordinates: fi, f*2, and R. The vector R 
connects the center of masses of the two HD molecules and is 
directed over the axis OZ, 9\ is the angle between fi and R, 62 
is the angle between F2 and R, ip2 is the torsional angle, ji,j2, 
and L are quantum angular momenta over the corresponding 
Jacobi coordinates fi, f*2, and R. 



II. METHOD: QUANTUM DYNAMICS 

Here we briefly present the close-coupling quantum- 
mechanical approach used in this study to calculate the 
cross sections and thermal rate coefficients of molecular 
hydrogen-hydrogen collision. The Schrodinger equation 
for the (12) + (34) collision in the center of a mass frame, 
where (12) and (34) are linear rigid rotors is [27l l28j: 
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where is the relative momentum operator, R is the 
relative position vector, M12 is the reduced mass of the 
pair M12 = (mi + m 2 )(m 3 + mi)/{mx + m 2 + m 3 + 
m i), A*i(2) are reduced masses of the targets: ^1(2) = 
m i(3) m 2(4)/("»i(3) + m 2(4)), h(2) arc tnc angles of orien- 
tation of rotors (12) and (34), respectively, E is the total 
ccnter-of-mass energy and V(fi ,?2,R) is the potential en- 
ergy surface for the four atomic system (12) + (34). The 
system is shown in Fig. [TJ Basically, the PESs of the 
H2-H2 and the HD-HD systems are the same. However, 
there is a small but important difference. To obtain the 
HD-HD PES from the existing H 2 -H 2 surface [H] one 
needs to appropriately shift the center of mass in the 
hydrogen molecules (H 2 ). The usual rigid rotor model 
[271 130| - |33| has also been applied in astrophysical calcu- 
lations of different atom and diatomic-molecule collisions 
or two diatomic-molecule collisions at low temperatures: 
T < 2000 K. 

The eigcnfunctions of the operators Lf 1{2) in Eq. ([T]) 
are simple spherical harmonics Yj imj (r). To solve Eq. ([1]) 
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where J is the total angular momentum quantum num- 
ber, M is its projection onto the space fixed z axis and 
the channel expansion functions are 
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with ji + j 2 = J12, ji2 + L = J, mi, 7712, mi2 and m the 
projections of ji, 32, 312 and L respectively. 

Substitution of Eq. © into ((T|) provides a set of cou- 
pled second order differential equations for the unknown 
radial functions U^ M (R) 



(- 



L(L+1) 
R 2 



+ K)u J a M {R) 



\dR 2 

2M isE / <€ M {h,h,R)\V{n,?2,R)\ 

a.' 

^J'ih^R) > U^, M {R)dhdf2dR, 



(4) 



where a = (jij 2 ji 2 i). We apply the hybrid modi- 
fied log-derivative- Airy propagator in the general purpose 
scattering program MOLSCAT [HI to solve the coupled 
Eqs. (U]). Additionally, we have tested other propaga- 
tor schemes included in MOLSCAT. Our calculations re- 
vealed that other propagators can also produce quite sta- 
ble results. 

Boothroyd et al. (BMKP) [29| constructed a global 
six-dimensional PES for two hydrogen molecules, espe- 
cially to represent the whole interaction region of the 
chemical reaction dynamics of the four-atomic system 
and to provide an accurate estimate of the van der Waals 
well. The ground state and a few excited-state ener- 
gies were calculated. In the six-dimensional configuration 
space of the H2-H2 system the conical intersection forms 
a complicated three-dimensional hyper surface. The new 
potential fits the van der Waals well to an accuracy of 
about 5% [H. In our calculation of the BMKP PES 
for H2+H2 the bond length was fixed at 1.449 a.u. or 
r£H 2 )=0.7668 A as in the Diep and Johnson (DJ) PES 
[3q . In the case of the HD+HD calculation the bond 
length of HD was adopted at r(HD) = 0.7631 A. 

The log-derivative matrix of the wave function is prop- 
agated to large i?-intermolecular distances, since all ex- 
perimentally observable quantum information about the 
collision is contained in the asymptotic behavior of func- 
tions U^ M (R — > 00). The numerical results are matched 
to the known asymptotic solution to derive the physical 
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scattering S'-matrix 
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where k aa i = \2M V2 (E + E a — E^)] 1 / 2 is the channel 
wave number, E a r a i^ arc rotational channel energies and 
E is the total energy in the (1234) system. The method 
was used for each partial wave until a converged cross 
section was obtained. 

Cross sections for rotational excitation and relaxation 
can be obtained directly from the ^-matrix. In partic- 
ular the cross sections for excitation from j\j2 — > j[j' 2 
summed over the hnal m\m'^ and averaged over the ini- 
tial TO1TO2 are given by 

tt(2J+ 1) 



x\6 aa> -S J U'i,j2,j'i2L';ji,32,ji2,L;E)\ 2 . (6) 



The kinetic energy is e = E- B 1 j 1 (j 1 + 1) - B 2 j 2 (32 + 1), 
where -61(2) are the rotation constants of rigid rotors (12) 
and (34) respectively. 

The relationship between the rate coefficient 
kj 1 j 2 ^.j' i j' 2 (T) and the corresponding cross section 
&j 1 j 2 ^j' 1 j!,( £ ) can be obtained through the following 
weighted average 
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where &g is Boltzmann constant, p is reduced mass of the 
molecule-molecule system and e s is the minimum kinetic 
energy for the levels j\ and j 2 to become accessible. 



III. RESULTS 

In this section our numerical results for rotational 
transitions in HD+HD collision and para/ para- hydrogen 
molecules are presented. We carry out state-to-state 
comparison between these two collisions for selected ro- 
tational transitions in the HD and H2 molecules. Specif- 
ically the following rotational energy transfer processes 
are considered: 



HD(ji) + HD(ja) -> HD(ji) 
Mh) + Mh) -> H 2 (i0 



-HD(j 2 ) 

W 2 ). 



(8) 
(9) 



At first look one might expect that the scattering out- 
puts of these two collisions © and © should be close 
to each other. This is because the PESs of H2-H 2 and 
HD-HD arc almost the same six-dimensional functions 
of the H4 four-atomic system coordinates. This fact 
follows from the general idea of the Born-Oppcnhcimer 



model [36j and simple theoretical atom-molecular consid- 
eration. Therefore, the two processes © and © should 
lead to similar results. At the same time the HD and H 2 
molecules have different rotational constants. This dif- 
ference is not dramatic: the rotational constant of H2 is 
B e (H.2) = 60.8 cm -1 , but the same parameter for HD is 
B e (HD) =44.7 cm" 1 . 

The HD+HD system has only four electrons. Fur- 
ther, the HD molecule consists of two nonidentical atoms 
which arc in a covalent bond. In covalent bonding the 
spins of the electrons are antiparallcl. The interaction of 
one of the nuclei, H + or D + , with its own electron leads 
to a quantum configuration where the spin of its electron 
is oriented antiparallcl to the spin of the nucleus, i.e. H + 
or D + . Thus, the spins of H + and D + are antiparallel. 
Because the spin of H + is I\ = —1/2 and the spin of D + 
is I2 = 1 the resulting spin of the HD molecule nucleus 
is /12 = 1/2. This value has been adopted in the current 
calculation, although there may be other possible values. 
The processes (JSj) and ((9]) arc collisions between two in- 
distinguishable diatomic molecules. This fact is taken 
into account in this computation. 



A. Comparison between HD+HD and H2+H2 
state-selected integral cross sections 

The precise HD-HD PES can be derived from the H 2 - 
H 2 surface by adjusting, i.e. shifting, the coordinates of 
the center of masses of the two H2 molecules to the center 
of masses of the HD molecules. Once the symmetry is 
broken in H2-H 2 by replacing the H atoms with the D 
atoms in the two H2 molecules we obtain the full HD- 
HD PES. The new potential will posses all parts of the 
HD-HD interaction. Therefore, it will be interesting to 
consider scatteing in two systems, which are not very 
different like HD+HD and H 2 +H 2 . 

In this work a large number of test calculations have 
been done to secure the convergence of the results with 
respect to all parameters that enter into the propagation 
of the Schrodingcr equation ([1]). This includes the inter- 
molecular distance R, the total angular momentum J of 
the four atomic system, the number of rotational levels 
Ni v i to be included in the close coupling expansion and 
some others, see the MOLSCAT manual [3J]. We reached 
convergence for the integral cross sections, a(Ekm), in 
all cases. However, it was particularly difficult to achieve 
convergence on the parameter R in both cases. For the 
applied BMKP PES we used R min = 2 A to R max = 50 
A. We also applied a few different propagators included 
in the MOLSCAT program. 

In a previous paper we presented a detailed description 
of convergence test for H2+H 2 collision [32|. The same 
numerical convergence has been achieved in this work. 
Namely, stable total cross sections have been obtained 
with respect to the number Ni v i of the rotational levels 
to be included in the basis set ((2]) of HD+HD, i. e., in 
each HD molecule the discrete integer quantum numbers 
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TABLE I: Rotational channel energies in the two hydrogen systems: a). HD + HD and b). para-H 2 + para-H.2 
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TABLE II: Comparison between different but "corresponding" (y a — vt) state-resolved cross-sections (A 2 ) in the HD + HD 
and para-K2 + para-H 2 collisions at ultracold T = 1.439 • 10 -8 K and very high T — 14390.0 K temperatures. 
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FIG. 2: (Color online) Integral cross sections for the HD+HD 
elastic scattering computed with the modified BMKP PES 
[2^ | together with the experimental and theoretical data from 
work 12011. 



ji and J2 run from to 4. As a result a maximum num- 
ber Ni v [(= 55) of rotational levels in HD+HD has been 
generated. With regard to the total quantum angular 
momentum J in the HD+HD system at T ~ 10~ 8 — 10~ 6 
K we needed to adopt just a few discrete values of this 
parameter, that is J ma x =3 or 4 was quite enough. How- 
ever, at larger collision energies E^mi such as, the ones 
we considered at T ~ 4000 K, one needs up to J max «60. 
In the case of E kin ~ 14000 K: J max rj 120. 

We present in Table I the rotational channel energies 
in the HD— HD and para-H2-para-H2 systems. This is 
a comparative table of the rotational spectra of these 
two systems. The first five columns from the left present 
HD— HD and the other five columns present the para- 
H2-para-H2 system; ji and ji are the quantum orbital 
momenta of the HD and H2 molecules, j'12 = j\ +32, 
with \j x - j 2 \ < J12 < ji +32, the index u a{b) is the cur- 
rent number of the degenerate rotational levels in HD and 
H2, respectively. The rotational energy levels are shown 
in cm -1 . The goal of this work is to investigate the ul- 
tracold regime in HD+HD, calculate its rotational energy 
transfer cross sections and thermal rate coefficients, and 
to carry out comparison with the corresponding (when 
v a = Vb, see Table I) rotational transitions in para-H2- 
para-H2 collision. 

There are two slightly different definitions [37j of 
the rotational cross sections in collisions between 
two identical diatomic molecules, for example, in 
HD(ji)+HD(j 2 H HD(ji)+HD(j 2 ). The cross section 
for the rotational transition jxj'2 —> is [27l . |34| : 

a ~ (1 + <5_ )1 j 2 )(l + Sf f ). However, in Ref. }38| the same 
cross section has been defined as a ~ (1 + SjifaSj/j/ ). It 
is seen that the two cross sections coincide when j\ + 1 j 2 
and j[ 7^ j 2 . However, for other combinations of the ro- 
tational quantum numbers, namely, when j\ = 32 and/or 
j\ = 79, the cross section calculated in accord with Rcfs. 
(271 [34J is two times larger than the cross section from 
Ref. 13811 . This has been taken into account in calculation 
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FIG. 3: (Color online) (a) Inelastic scattering integral cross 
sections for HD(0) + HD(1) -> HD(0) + HD(0) and H 2 (0) + 
H 2 (2) H 2 (0) + H 2 (0). (b) The same for HD(1) + HD(2) -)■ 
HD(1) + HD(1) and H 2 (2) + H 2 (4) -)• H 2 (2) + H 2 (2). 



with the MOLSCAT program [3J|, i.e., for the integral 
cross sections cr(j' 1 j' 2 ;jij2,e) from Eq. ([6|) the following 
pre-factor [(1 + 6j t j 2 ) (1 + has been adopted. 

First, let us turn to HD+HD elastic scattering. In 
Fig. [5] we show results of this work computed with the 
modified BMKP PES [2!| together with the correspond- 
ing data (experiment/theory) from relatively old papers 
[201 ] . As can be seen all these cross sections are in a 
satisfcatory agreement with each other. This test calcu- 
lation reveals the reliability of the modified BMKP PES, 
the computer program and the numerical convergence. 
One can see that at low and very low energies the gen- 
eral forms of the cross sections are rather close to each 
other with the exception of a shape resonance and small 
oscillations in the cross sections in Fig. 2. 

In Table II few selected state-to-state cross sections for 
the HD+HD and H 2 +H2 collisions are presented. Here 
we compare results for few specific rotational excitation 
and dc-excitation integral cross sections at only two val- 
ues of kinetic energy, namely at T\ ~ 10~ 8 K, and at 
very high collision energy, i.e. T 2 ~ 14, 000 K. At low 
temperature the HD+HD cross section could be larger 
by three to four orders of magnitude, whereas at high 
temperature the two cross sections are of the same order. 

In Figs. [31 Q] and [5] we present a few state-selected 
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FIG. 4: (Color online) (a) Inelastic scattering integral cross 
sections for HD(1) + HD(2) -> HD(0) + HD(1) and H 2 (2) + 
H 2 (4) -J- H 2 (0) + H 2 (2). (b) The same for HD(2) + HD(2) -> 
HD(1) + HD(2) and H 2 (4) + H 2 (4) -> H 2 (2) + H 2 (4). 



FIG. 5: (Color online) (a) Inelastic scattering integral cross 
sections for HD(2) + HD(2) -> HD(1) + HD(1) and H 2 (4) + 
H 2 (4) ->■ H 2 (2) + H 2 (2). (b) The same for HD(2) + HD(2) -)■ 
HD(0) + HD(2) and H 2 (4) + H 2 (4) -)• H 2 (0) + H 2 (4). 



rotational transition cross sections in collisions ([8j and 
©. It is useful to see the corresponding cross sections 
together on a single plot, that is when v a = v\,. For 
example, in Fig. [3] (a) we show rotational transition de- 
excitation cross sections from the first excited states of 
HD and H2 molecules, i.e. we consider HD(0)+HD(1) — >• 
HD(0) + HD(0) and H 2 (0) + H 2 (2) ^H 2 (0) + H 2 (0) for a 
wide range of kinetic energies: from 1 K to up to 4000 K. 
In the case of H 2 +H 2 , we carry out computations with 
two different PESs, e.g., with the BMKP PES [H| and 
with the Diep- Johnson (DJ) H 2 -H 2 PES from Ref. 
The last one was formulated for fixed equilibrium dis- 
tances between the hydrogen atoms in each H 2 molecule. 
In Fig. [3] (b) we show cross sections for some other de- 
excitation processes in the HD-HD system. 

Further results for the integral cross section are shown 
in Figs. 0] and [SJ The cross sections for following pro- 
cesses are presented in Fig. [4] (a): HD(1) + HD(2) -)• 
HD(0) + HD(1) and H 2 (2) + H 2 (4) H 2 (0) + H 2 (2). 
The cross sections of HD(2) + HD(2) -> HD(1) + HD(2) 
and H 2 (4) + H 2 (4) -> H 2 (2) + H 2 (4) are shown in Fig. 
H (b). It is seen in Fig. |J (a) that the HD+HD de- 
excitation cross sections could be larger than the H 2 +H 2 
de-excitation cross sections by four orders of magnitudes. 
For the cross sections in Fig. H](b) the difference is about 



a factor of 10 3 . In Fig. [5] (a) we show results for rota- 
tional transitions: HD(2) + HD(2) -> HD(1) + HD(1) and 
H a (4) + H a (4) -> H 2 (2) + H 2 (2). In Fig. E (b), we show 
the same for the following transitions. HD(2) + HD(2) — > 
HD(0) + HD(2) and H 2 (4) + H 2 (4) -»• H 2 (0) + H 2 (4). 
Again, in both cases the HD+HD cross sections could be 
larger than the H 2 +H 2 at 1 K by four orders of mag- 
nitudes. We see in Figs. 3, 4, and 5 that, for the de- 
excitation processes for H 2 +H 2 , the BMKP and the DJ 
PESs provide similar results with the exceptions shown 
in Fig. 2. 

In Ref. [39l | it was shown that for a specific excita- 
tion rotational transition in the H 2 +H 2 inelastic scatter- 
ing, e.g., in H 2 (0) + H 2 (0) -> H 2 (6) + H 2 (2), the BMKP 
PES provides an incorrect cross section when compared 
to the DJ potential. The comparison was also carried 
out with available experimental data |4(|. Neverthe- 
less, the BMKP PES has been applied to the important 
(astrophysical) o-/p-H 2 +HD inelastic scattering problem 
[HI, Ha). This is why we applied the BMKP PES to the 
HD+HD scattering problem. 

In conclusion, based on the Born-Oppcnhcimer model 
treatment, the HD+HD and H 2 +H 2 have similar PESs. 
However, in the case of the HD-HD system the origi- 
nal H 2 -H 2 PES is adopted and the two center of masses 
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FIG. 6: (Color online) Thermal rate coefficients for the inelas- 
tic scattering processes: HD(0) + HD(1) -> HD(0) + HD(0), 
HD(1) + HD(1) -> HD(0) + HD(1), and HD(1) + HD(1) ->■ 
HD(0) + HD(0), at ultracold temperatures. 



FIG. 8: (Color online) The same as in Fig. |S]for following pro- 
cesses: HD(1)+HD(2) -> HD(1) + HD(1), HD(1) +HD(2) -s- 
HD(0) + HD(2), HD(1) + HD(2) -> HD(0) + HD(1), and 
HD(1) + HD(2) -> HD(0) + HD(0). 



of both H2 molecules is just slightly shifted to the ap- 
propriate positions of the HD molecule center of masses. 
After this procedure we obtain the full space, i.e. global 
HD-HD PES. Our computations with this modified PES 
revealed a very strong isotopic effect in the HD+HD and 
H2+H2 collisions at low energies. 



B. HD+HD rotational state-selected thermal rate 
coefficients at ultracold temperatures 

We show in Figs. El El El an d M the thermal rate co- 
efficients in the inelastic HD-f-HD collision at very low 
temperatures from ~5x 10~ 7 K to ~ 2 x 1CT 5 K. These 
results were obtained from corresponding state-resolved 
integral cross sections ^j 1 j 2 ^j' 1 j 2 ( e ) with the use of ex- 
pression ([7]) . Only de-excitation thermal rates have been 
computed, because at such a low temperature the exci- 
tation thermal rates are extremely small. The rates have 



-HD(O) + HD(2) HD(O) + HD(1) 
1 HD(O) I HD(2) > HD(1) I HD(1) 
■ HD(O) + HD(2) -+ HD(O) + HD(O) 



Temperature T (K) 

FIG. 7: (Color online) The same as in Fig. [6]for following pro- 
cesses: HD(0) + HD(2) -> HD(0) + HD(1), HD(0) + HD(2) -> 
HD(1) + HD(1), and HD(0) + HD(2) -> HD(0) + HD(0). 



been computed for different initial rotational states of 
the HD molecules. The figure captions include the infor- 
mation about the specific state-selected rotational tran- 
sitions in both HD molecules before and after collision. 

In Fig. [SI we show results for the rates kj 1 j 2 ^j> i j 2 (T) 
for the inelastic scattering processes from the first three 
rotational states. In Fig. [7] we show the rates for tran- 
sition from the excited states. These transitions corre- 
spond to different de-excitations from lower initial ro- 
tational states in the HD— HD system. In Fig. [6] the 
initial state for the solid line corresponds to the first ex- 
cited rotational state in HD+HD. From Table I one can 
see that for this state of the system: ^ a =2. The two 
other lines in Fig. [6] correspond to the initial state with 
v a =A : i.e. £^ = 178.8 cm -1 . The rates for the initial 
state HD(1)+HD(1) in Fig. [Jo] are smaller than those for 
the initial state HD(0)+HD(1), although all three rates 
have similar dependence on temperature. 

In Fig. [7| we show the rates for some de-excitation 
transitions from the initial state HD(0)+HD(2), which is 
the third rotational excited state in HD+HD: ^a=3. Al- 
though we again observe that the behavior of the thermal 
rates k 



(T) is quite identical, their values signif- 
icantly differ from each other, specifically up to two or- 
ders of magnitude. Further, the thermal rate coefficients 
kj 1 j 2 -^j' j 2 (T) from the higher excited rotational states 
HD(1)+HD(2) and HD(2)+HD(2) are presented in Fig. M 
and respectively. Specifically, these rates are from the 
energy levels 357.6 cm -1 and 536.4 cm -1 corresponding 
to the following two indices: u a =b and ^ a =6, respectively 
(see Table I). In these calculations we needed a fairly ex- 
tended number of basis functions in the expansion @ 
for convergence. However, this is quite understandable 
because of the large energy of the initial state. For exam- 
ple, in the calculation of rotational transitions from that 
level with the following rotational indices j\ = 2, ji = 2, 
i.e. v a = 6 and =536.4 cm -1 (Table I) all lower lying 
rotational levels have to be included in the computation. 
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FIG. 9: (Color online) The same as in Fig. [S] for follow- 
ing processes: HD(2) + HD(2) HD(1) + HD(2), HD(2) + 
HD(2) -> HD(1)+HD(1), HD(2)+HD(2) HD(0) + HD(2), 
HD(2) + HD(2) HD(0) + HD(1), and HD(2) + HD(2) -> 
HD(0) + HD(0). 



It is known that a quantum-mechanical transition 
probability P a ^p between any two quantum states, for 
example, the initial a — (jij'2) and final f3 ~ (j'ij'2), is 
inversely proportional to the energy gap Ae a p between 
these two states. In turn, the cross sections a a p(E) 
and corresponding thermal rate coefficients fc Q ^(T) are 
directly proportional to the quantum probabilities. Thus 



(10) 



However, a quite unexpected result relating the rates in 
the following reactions is seen in the present study: 



HD(0) 
HD(0) 



HD(2) 
HD(2) 



HD(0) 
HD(1) 



HD(1), 
HD(1). 



(11) 
(12) 



It is seen from Fig. [7J that process (fTTj) has much 
larger thermal rate coefficients than process (fT2"j). The 
difference between these rates is about an order of 
magnitude. At the same time the energy difference 
between the HD molecule rotational states in (fTTj) is 
Aeo2-oi=178.8 cm -1 , which is larger than in (fT2"j) . for 
which Aeo2-n=89.4 cm -1 . In accord with the relation- 
ship ([TO)) one could expect that the process (fTTj) would 
have lower values of the thermal rates than the process 
(fT2"j) . but it does not. This happens because in ([T2"j) 
both HD molecules simultaneously change their internal 
states, i.e. rotational quantum numbers. Probably, this 
is the reason that the process fT2j) is much slower than 
(fTTj) . This result is somewhat similar to older computa- 
tional data on the H2+H2 collision (32l.|43j. where the au- 
thors found that the excitation process H2 (0) + H2 (0) — >• 
H2(4) +H2(2) has larger cross sections at larger collision 
energies than the process H 2 (0) + H 2 (0) -> H 2 (4)+H 2 (0). 



IV. CONCLUSIONS AND FUTURE WORK 

Currently theoretical and experimental research in the 
field of the molecular Bose-Einstein condensates at ul- 
tracold temperatures is increasingly gaining momentum 
0, HH, |45j|. For example, in the recent work [4|| the 
authors develop a promising approach for laser cooling 
of diatomic polar molecules. The method should allow 
the production of large samples of molecules at ultracold 
temperatures. Only a few of the possible practical and 
technological applications where new results of this re- 
search could be used have been briefly outlined in the 
Introduction. Researchers in this new field of atomic, 
molecular and optical physics have had tremendous suc- 
cess within last two decades. It is useful to have ex- 
act, high quality full space PESs for such polar molecule 
interactions. In turn the HD+HD system could be a 
prototype collision between two polar molecules with a 
high quality full space four- atomic PES. In this work 
we performed a detailed quantum-mechanical study of 
the state-resolved rotational excitation/de-excitation col- 
lisions between hydrogen molecules. The HD + HD — > 
HD+HD and H2+H2 — > H2+H2 collisions have been con- 
sidered and their rotational state-selected integral cross 
sections have been computed for a wide range of temper- 
atures, i.e. from ultracold T ~ 10~ 8 K to up to T ~ 
14000 K. We have demonstrated that a small change in 
the H 2 -H 2 PES to adjust for HD-HD can lead to sub- 
stantial differences in the scattering outputs, i.e. in the 
integral state resolved cross sections. This calculation 
was carried out within a single H4 PES from Ref . [29[ . 

Further, in connection with the problems of coherent 
control of the atomic and molecular interactions the au- 
thors of Ref. [Ill performed a numerical investigation 
of the quantum entanglement for the case of a nonre- 
active ultracold collision between two indistinguishable 
(H2+H2) molecules. Similarly, the initial state quan- 
tum entanglement coupling has been considered for the 
case of ultracold collision between identical two-atomic 
polar molecules [13] . It would also be useful to men- 
tion here that universal relations for strongly correlated 
fermions have been derived recently [48|, |4i| . Because the 
HD molecules are fermions with a well known interaction 
potential, they could be useful, for example, in direct nu- 
merical verification of these universal relations. 

The authors of Ref. [5(| formulated a time- 
independent quantum-mechanical formalism to describe 
the dynamics of molecules with permanent electric dipolc 
moments in a two-dimensional confined geometry such as 
in a one-dimensional optical lattice. It would be useful 
in future investigation to adopt these techniques and ap- 
ply them to a fermionic molecule system such as HD+HD 
with a well-known potential [29j.[5l|. It would be interest- 
ing to see differences in the quantum dynamics between 
the state-resolved HD+HD rotational thermal rate co- 
efficients of the current work and possible new thermal 
rates for the same system but when embedded in a one- 
dimensional optical lattice or a microwave trap. 
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